TWO VARIABLE DEFORMATIONS OF THE CHEBYSHEV 

MEASURE 

JEFFREY S. GERONIMO AND PLAMEN ILIEV 



Abstract. We construct one and two parameter deformations of the two di- 
mensional Chebyshev polynomials with simple recurrence coefficients, follow- 
ing the algorithm in [3]. Using inverse scattering techniques, we compute the 
corresponding orthogonality measures. 



1. Introduction 

The general theory of bivariate orthogonal polynomials goes back to the work of 
Jackson [T2]. Special examples have arisen in studies related to symmetric groups 
[HEIl [22]), as extensions of one variable polynomials and as eigenfunctions 

of partial differential equations [IH1 HH H31 [21] ■ An updated account of the theory 
can be found in the books [5J [25] . 

A major difficulty encountered in the theory of orthogonal polynomials of more 
than one variable is which monomial ordering to use. Starting with [12] . the pre- 
ferred ordering is the total degree ordering and for polynomials with the same total 
degree the ordering is lexicographical that is 

(k, I) <totaldeg (kl,h) 

if 

k + I < ki + h or (k + I = h + h (k,l) <i ox (ki,h)), 
where lex means the lexicographical ordering (see below). In other words, we apply 
the Gram-Schmidt process to the polynomials ordered as follows {1, y, x, y 2 , xy, x 2 , . . 
If we denote by P n (x,y) the (n + 1) dimensional vector with components the or- 
thonormal polynomials of total degree n, then the multiplications by x and y are 
given by three term recurrence relations 

xP n = A x , n P n+ i + B x>n P n + A^ n _ 1 P„_ 1 (1.1) 

yP n = A y , n P n+1 + B v . n P n + 4,n-iP»-ii (1-2) 

where A XjTl , A y _ n are (n+1) x (n + 2) matrices such that rank(A :Ei „) = rank(Aj /i „) = 
n + 1 and B x n , B Vtn are symmetric (n + 1) x (n + 1) matrices. Notice that the 
Jacobi matrices corresponding to multiplications by x and y will commute which 
amounts to certain commutativity relations between the matrices defined above, 
see [SJ for more details. 

Inspired by the progress made in [10) in connection with Fejer-Riesz factoriza- 
tions and autoregressive filters in two variables, an alternative way to approach 
two dimensional orthogonal polynomials was proposed in [3] by relating them to 
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the theory of matrix valued orthogonal polynomials. This can be accomplished by 
using the lexicographical ordering, 

(k,l) <ie X (ki,h) <^> k < ki or (k = ki and I < li), 

or the reverse lexicographical ordering 

(k,l) <revlcx (h,h) ^ (l,k) <lex (h,h), 

to arrange the monomials. This naturally connects the theory of bivariate orthog- 
onal polynomials to doubly Hankel matrices. In particular, this led to an alternate 
parametrization for positive doubly Hankel matrices and an algorithm for construct- 
ing their inverses. 

In the present paper we apply the algorithm from [3] to construct deformations 
of the two dimensional Chebyshev polynomials with relatively simple recurrence 
coefficients. These polynomials can be thought of as two dimensional analogs of the 
Bernstein-Szego polynomials. Using the connection of the above theory with matrix 
valued orthogonal polynomials and the inverse scattering techniques developed in 
[9j we show how one can obtain the orthogonality measure from the recurrence 
coefficients. 

The paper is organized as follows. In the next section we recall the main ingre- 
dients needed for the construction. In particular, we review the parametrization 
of the doubly Hankel matrices in [3] and the algorithm which gives the recurrence 
coefficients. We also sketch the inverse scattering techniques in the matrix case 
and the connection with the Darboux transformation. In Sections 3 and 4, respec- 
tively we present one parameter and two parameter deformations of the recurrence 
coefficients for the two dimensional Chebyshev polynomials and we derive the or- 
thogonality measure. 



2. Bivariate Orthogonal Polynomials 
Let /j,(x,y) be a Borel measure supported on R 2 , such that 



f(x,y)dfi(x,y) < oo 



for every polynomial f(x, y). For every nonnegative integer m we order the mono- 
mials x'y 3 , < i < n, < j < m lexicographically i.e. 



■■■,y m ,x,xy, ...,xy m ,x 2 , 



}• 



We construct the moment matrix 



H Hi 

Hi H 2 



H n H n+ i 



(2.1) 
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where each Hi is an (to + 1) x (to + 1) matrix of the form 
hn hi2 



Hi 



'-■2.2 m, 



0, . ..,2n. 



(2.2) 



Thus H niJn is a block Hankel matrix where each block is a Hankcl matrix so it has 
a doubly Hankcl structure. If the reverse lexicographical ordering is used in place 
of the lexicographical ordering we obtain another moment matrix H n ^ m where the 
roles of n and m are interchanged. We will assume that H n ^ m is positive definite for 
all n and to. For every nonnegativc integer to we apply the Gram-Schmidt process 
to the basis of monomials ordered as above and define the orthonormal polynomials 
Pn m( s i v)i < n, < to, < ^ < to, by the equations, 



/ 
/ 



Pn m^y^nix, y) = 0, < i < n and < j < to or z — n and < j < I, 



P l n. m Pn, m Mx,y) 



1, 



and 



P l n , m {x,y) = k^ l x n y l + Yl 

(*,j')<lex(n,0 



K n,m,l X V ■ 



(2.3) 



(2.4) 



With the convention fc™' TO ; > 0, the above equations uniquely specify p l n m . Poly- 
nomials orthonormal with respect to dfi(x,y) but using the reverse lexicographical 
ordering will be denoted by p l nm - They are uniquely determined by the above 
relations with the roles of n and to interchanged. Set, 



f n, m 

p 1 

i n.rn 



V 



(2.5) 



The P n ,m uiay be obtained in an alternate manner as follows. We associate an 
(to + 1) x (to + 1) matrix valued measure dM m+1 (x) by taking 

dM m+1 (x) = [ [l,y, . . . , y m ]*d^(x, y)[l,y, y m ], (2.6) 

where the above integral is with respect to y. Let us denote by {P n ,m}n°=o the 
sequence of (to + 1) x (to + 1) matrix valued polynomials satisfying 



/ 

JR 



P n ,m{x) = K™ n x n + lower order terms, 
P^ m (x)dM m+1 (P Km (x)Y = S k , m I m+1 , 



(2.7) 
(2.8) 



with K™ n a lower triangular matrix with strictly positive diagonal entries. The 
above conditions uniquely specify these left matrix valued orthogonal polynomials 
and it follows that 

F„, m (s, y) = P n , m {x)[l,y, y m Y. (2.9) 
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From equation (12. 6| we see that M m+1 is a Hankel matrix and we have, 

Theorem 2.1. Let M rn+1 (x) be a matrix measure supported on some interval 
I C M then M m+1 (x) is a Hankel measure if and only if there exists a bounded 
positive linear functional S^m : C(I) x {1, . . . , y 2m } — > K with (J 7 x k dM m+1 )ij = 

Proof. If we begin with M m+1 (x) set H, = JjX l dM m+1 and let H n>m be the 
(n + l)(m + 1) x (n + l)(m + 1) matrix given by equation (|2.ip . Then H njlri 
is a positive definite doubly Hankel matrix which implies via the above formula 
connecting M m+1 and U^m and the density of the polynomials in C(I) that 
is a bounded positive linear functional which is uniquely defined. The converse 
follows from the relation between ?2m an d M m+1 given above. □ 

The following Theorem was proved in [3], 

Theorem 2.2. Given {P rl , TO } and {P raim } the following recurrence formulas hold, 

(2.11) 
(2.12) 
(2.13) 



r t 


D _ Tp 


1 ^71,771^71— 1,771 7 


J 1 1 


n.rn — y^n,m- 


-1 H~ Jn^jr^n—l^m H~ ^n.rn^n 




— J TP 4 


" ,771-^^,^—1 




^n,m — 


{x^n — l,m) -^71,771)) 










/i - 


(yPn, 771 — 1) -^r^TTi) 




7 2 - 

^ n,m 


— (yIP n _ m _ 1 , IPn — l,m) 






— (y 1^71 ,?7i - 1 > IPn - 1 , 771 - 1 ) 




r — 

J- n,m — 


(-^71,771 — 1) -Pn,77l) 




, rn — 


(-^71,771— 1) ^ti— l,m) 




, ? n — 


(-Pn,771 ) ^71, 771 ) 



where for every two vector valued polynomials /, 5 we sef 

(/,<?) = / f{x,y)g t {x,y)d^x,y). (2.14) 



Similar formulas hold for ¥ n Tn (x, y). 

Equation (|2. 10|) comes from the fact that the vector polynomials can be related 
to matrix orthogonal polynomials as discussed above. Given the coefficients in this 
equation it is possible to compute the vector polynomials in a strip of size m + 1 
in the n direction. Its tilde counterpart obtained from the reverse lexicographical 
ordering allows one to compute the tilde vector polynomials along a strip of size 
71 + 1 in the m direction. Equations (|2.1ip and (|2.12[) allow one to compute by 
changing n and m simultaneously. That is given the polynomials at the (n, m — 1) 
level and the (n — l,m) level and the coefficients in ([2. lip . (|2.12|) and (|2. 13|) it is 
possible to compute V 7hm and P n ,m- It should be noted that in going from (n— 1, m), 
and (n, m — 1) to (n, m) there are generically 4 new moments. Examination of the 
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sizes of the coefficients in (|2.1ip and (|2.12[) shows that most of the coefficients in 
these matrices must be computable from the coefficients given at previous levels. 
Indeed in [3] it was shown that this is the case and that for n > 1 and m > 1 only 
(Xn,m)m,n, (^, m )m,n> ( J n, m )m,m and (J^ m ) m ,m+i need be specified in order to 
compute P n ,m- For n > 1 and m > 1 we denote 

(3^-n,m)m,n = s 2n — 1,2m— 1 m )m,n = s 2n — 1,2m (2.15) 

(.^n,m)m : m — ^277,2777—1 (^n.m/^.^i+l 3277,2777- (^'1^0 

For n — and m > or n > 1 and m = we introduce sq m and s n! o which 
are parameters in the one dimensional recurrence formulas associated with the line 
(n, 0) or (0, m) respectively. Thus we find 

Theorem 2.3. [3] Given parameters sq.o, ■ ■ ■ , S2n.2m S R, we construct 

• scalars A^ifi, B^o, i = 0, . . . ,n — 1, and j4o,j+i, -Boj, j = 0, . . . ,m — 1/ 

• j x i matrices DQ j and J?,, i = 1, . . . ,n, j = 1, . . . , m; 



3 x (j + 1) matrices J\ ■ /or £ = 1, . . . , n } j = 1, 2 . 



, m. 



s 2i ,2i > and | < 1, (2.17) 
i/ien iftere exists a positive linear functional 5" smc/i £/ia£ 

S'(Pi,m,Pj,in) = JijUl and J(Pn,„ Pnj) = d^n+L (2.18) 

TTie conditions (|2. 1T[) are aZso necessary. 

With the above we are able to begin with the parameters Sij and try to compute 
two variable orthogonality measures associated with these parameters. The con- 
nection with matrix orthogonal polynomials will be especially useful and we review 
some of the relevant relations. The theory of matrix orthogonal polynomials has 
an extensive literature (TJ [6l [9l [23l [24l [26] and we begin with the equation (|2.10l) 
and use (|2.9| to find, 

%Pn,m ^77+1, 777-^77 + 1, 777 ~T" -^77,777-^77,777 ^77 ,777 -^77 — 1 ,777 • (^'^^) 

If we set -P-i.7 n = and Po. m = Im+i then the above equation uniquely gives 
{Pn,m}- We now turn the problem around and begin with {A n/m \ ( ^' =1 and {-B„. m }^L 
where Aj iTn Bi_i. m are (m + 1) x (m + 1) real valued matrices with B n/m = B l n m 
and An,?™ lower triangular with strictly positive diagonal entries. We will suppose 
that 

lim A n m = -I m +i and lim B n m = 0, (2.20) 

77 — >oo Z n — y oo 

and define Ao :m = I m +i- In this case, there is a unique matrix measure M such 
that 

J P n ,m{x)dM{x)P ktm {x) t =S ntk I m+1 . (2.21) 

Following [9] we introduce the matrix valued function 

^ n,m{.z) — Pn.m, 2zA nrn P n _i m , (2.22) 



where z — x — y 'x 2 — 1. The branch of the square root is chosen so that z — > as 
x — > +00. Using (|2.19[) and ()2 . 22[) one can deduce that 

1 _ 1 _ t (2.23) 

^^^77, 777 ^ 77— 1,777 "(" ^^"'l,™ 4yl„ , m A n n% ) Z — 2B n _l m ^ -P77— 1,777' 
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Note that if A n>m = ^I m +i and B„_i >m = for all n> uq then 

*; m (2) = z n ^ n ,n{z) = *; 0>ro (^) for all n > n . 
We also introduce two scattering solutions of (|2. 19[) P^ m {z) that satisfy 

lim ||^"P± (z)-/ m+1 ||=0 



where the matrix norm is the Hilbert- Schmidt norm i.e \ \B\ \ = (ti(BB^)) 1 ^ 2 . If we 
assume 

oo 

^n(||l-4^ n , mJ 4^ m || + ||S n _ 1)m ||) <oo, (2.24) 

71=1 

then a slight modification of the techniques in [9] gives that for n > 0, P^~ m exists, 
is continuous for \z\ < 1 and analytic for \z\ < 1 while P^m exists, is continuous 
for \z\ > 1 and analytic for \z\ > 1. Since the coefficients in (12.19)) are real valued 
we see from the asymptotic conditions satisfied by P^m that for \z\ — 1, P,7m( z ) = 
P+ TO (l/z). For two solutions Y and X of (|2.19[) we define 

W m [X,Y] = Xl m (x)A n+hm Y n+hm (x) - X^ m (z)^ +1)ro Y n>ro (z), (2.25) 

which by standard computations is independent of n. With the use of W m it follows 
that if ([HID holds then for all \z\ = 1, z ^ ±1 

= ^7^ P "V (*)/-(*) - (2-26) 



where 
and 



/+(z) = P+ 1 , m (z) t |z|<l, (2.27) 
/+(*) = ^ I™ *U(*)> N < 1- (2-28) 



From the relation between P^m we find, 

/+(*) = /-(l/*) \z\ = l (2.29) 
Manipulations similar to those that lead to (12.25)) yield 

1 (2-30) 
= -(z-l/z)I m+1 , \z\ = l, z?±l, 

which implies that P + is nonsingular for \z\ — l,z ^= ±1. This follows since if 
there is a vector a and a zoil^ol — 1j z o 7^ il such that P+ m (zo)a = then 
(P+ m (z ) a ) t = which cannot happen by the equation above. This implies that / + 
is nonsingular for \z\ — 1, z ^ ±1. With n = —1 in (|2 . 26)) we find for \z\ = 1, ±1 
the useful equation 

/+/- = /-/+• (2-31) 
Another useful relation between / + and P+ m is [HI formula 7.3], 

P+ m (z) = ^ EnfM dM {y)f + {z) \z\ < 1 (2.32) 
J-oo x-y 

where the relation between x and z is given above. 
We will now make the following assumptions 

A n +i,m = 2^m+i) Pn,m = 0, Vn > n , (2.33) 
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and 

det(z/+(z)) ± 0, \z\ < 1, (2.34) 

which will simplify the presentation below and are sufficiently general for our exam- 
ples. The more general case will be taken up later (see however [9] and [21] )• With 
assumption ([2~33]) we see from (|2~28| and (|2~23f that 2zf + (z) = V* , m {z) which 
will be a matrix polynomial of degree 2uq if A n0:7n ^ |im+i an d of degree 2uq — 1 
if A no , m = i/ m+ i and _B„ _i. m 7^ and for n > 0, -P^ m are matrix polynomials in 
z^ 1 respectively. Furthermore (]2 . 26[) and contour integration show that the matrix 
orthogonal polynomials {P n ,m\ satisfy equation (|2.21[) with 

dM(x) = a m (x)dx = — \J\ - x 2 (f+(z) f f+(z))- l dx, z = e w . (2.35) 
2n 

What we would like to calculate next is what happens when we add a mass point 
to the system which can be accomplished by multiplying ^n .m by a zero. To this 
end write, 

*„ +i,mO) = (z /z - l)* n0im (z), (2.36) 

where z is real and \zq\ < 1. With the choice of $ no+ i tTn the computation of 
the orthogonality measure as well as the matrix orthogonal polynomials {P n ,m(x)} 
with -Po,m = Im+i can be carried out in the following manner, 

Theorem 2.4. Given {A n ^ m } and {B n m } satisfying (|2 . 33[) and ^ no ,m satisfy- 
ing (|2.34p . Let i& no +x tm be given by (|2.36|) . Then {P n .m(x)} are a set of matrix 
orthogonal polynomials satisfying 



P n , m (x)dM(x)P k , m (x) t dx = I m +iS n ,k n,k>0 (2.37) 



where 



with 



dM(x) = {P^> \ <X< d ) , ( 2 - 38 ) 
[r m d(x - x ), x f [-1, 1\ 

f m (i) = ^ (2-39) 
2z 

and 

4z 

The polynomials {P n , m } have the Uvarov-Christoffel representation, 

-Pn,m(^) — dn,m(Qn — l,m(xQ^ A n ^ m P n ^ m (x^ Qn.m A n m Pn — l,m (^)) 3 (2-41) 

where 

/>oo 

(x ) = / P n , m {x)dM{x). (2.42) 



The constants d n ^ m are chosen so that P n , m is orthonormal for n > 0, which is 
equivalent to 

(di^d^)- 1 = ^l^/_(z )- 1 p-_ l m (z ) t A„,„ l p- m (z )/i(z )- 1 . (2.43) 
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Proof. From the definition of Q n ,m(x), dM(x) and a m we find for n > and 
< i < n that 



P n , m {x){x - x ) l dM(x) = 0. (2.44) 

— OO 

Equations f2~32"]) . (pOTj) and fUgJ show that 

Q n (x ) = P+ m (z )/|(z )- 1 /2z + P n , m ( 2 ;o)f m = P- m (z )f_ {z^ 1 /2z Ql (2.45) 

where we have used the definition of f m and the analytic properties of P^m- We use 
the above equation to define Q_\ m . The substitution of this formula into (|2.44[) 
with i — yields 

/OO 
P n , m (x)dM(x) 
-oo 

= (2z ) 2 ^~( Z °^ t ^n-X,m{ z ^) t ^n,m,Pn,m( z ^) ~ ^n,ro( z o)*^n, m ^_i, m (2o))/+(2 ; o) 1 
+ ^{f-izayHPn-l^ofAn^Pn^ixo) - P~ m (z )* A^ m P„_l >m (x ))f m . 

(2.46) 

If in the first term on the left hand side we use (|2.25J) and let n tend to oo while in 
the second term we use (|2.25|) and set n = 0we find from the definition of f m and 
/+ that the above integral is equal to for n > 0. To show that d n , m is given as in 
equation (|2.43|) we begin with the recurrence formula satisfied by Q n ,m(xo) 

XoQn,m — ,mQ n+1 ,m ~\~ B n m Q n rn -\- -A. n m Q n —i m7 (2.47) 

which follows from (|2.42[) . Thus routine manipulations give 

d n}nPn, m (x) = (X- .T ) ^ Ql : ,n( X o)Pi;rn(x) + Q~l, m {x ), 
i=0 

so that for n > 

P n , m (x)dM{x)Pl m (x) 

-ir-d n , m y2Qt {z ) / P iim (x)dM{x)Pl m (x) 
-z d n ^ m Q n _ 1 (zo)A n ^ m Q ntm (zo)d n . 

ZZq 



(2.48) 



The last equality follows by eliminating P n ^ m using its definition. □ 

Remark 2.5. If we want to fix d n ^ m so that the highest coefficient of P n ,m(x) 
is a lower triangular matrix, we need to choose in (|2.43|) the unique solution for 
which dn^mQn-i^ixo)^ is a lower triangular matrix. Notice that for n > no we 
have P~ m (z ) = z„ n I m+ i, A n+ i <m = \l m +i which leads to d n+1 . m = Az^ +2 f-(z ). 
Thus for n > no, formula (|2.41l) reduces to P n , m — zoPn,m ~ Pn—i,m- This implies 
A n +i, m = \l m +i and B n . m = for n > n . 



Remark 2.6. Theorem l2.4l can be easily explained in terms of the Darboux trans- 
formation. Indeed, let SL m be the second order difference operator corresponding 
to the second order difference equation (|2.19p for the polynomials P n ,m{%), he. 
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(&m(f))n = A n+ i <m f n +i + B n , m fn + A^ TO / n _i. We can think of L m as a second 
order difference operator acting on matrix valued functions (the size of the matrices 
is (to + 1) x (m + 1)) of a discrete variable n. Similarly, we denote by L m the op- 
erator corresponding to the polynomials P n ,m{x). Then using (|2.47|) one can check 
that operator L m — x Id can be factored as 

£ m - x Id = TQ, (2.49) 

where Q is the backward difference operator in (|2.41j) , i.e. 

n— l,m (X ) t < ro /n-l) (2.50) 

and T is the forward difference operator 

•P(/)n — Wn,m) (/n+l^n+l,m — /n^n,m) * 

The operator L m is obtained from L m by exchanging the factors in (|2.49[) . i.e. 

L m - x Id = Q3>. (2.51) 

The above discussion can be summarized as, 

Theorem 2.7. Suppose that the coefficients in L m satisfy (|2.33[) and zf+(z) sat- 
isfies (|2.34p . If L rn is related to L m by the Darboux transformation (I2.49P and 
(I2.5ip and Q is a backward difference operator of the form (|2.50[) where Q n ,m{xo) 
are given by (|2.45[) i/ien t/ie polynomials P n _ m are orthogonal with respect to M 
which is given by equations ([2~37|) - (j^40|) . 

We can now use the above results to help solve the bivariate problem via the 
parametric moment problem. Equation (|2.6p and Theorem 1 2. II show that dM rn+1 is 
a Hankel matrix with entries dfij(x) — J y°dli{x, y) with the integration being over 
y. The above construction allows us to compute Hj(x) for every j. We can think of 
Hj(x) as one dimensional moments (in the variable y), depending on a parameter x. 
Let us consider the corresponding polynomials q^ l (y). The three term recurrence 
formula takes the form 

with coefficients depending also on the parameter x. Next we use one dimensional 
theory and compute dfi(x,y) by introducing the function 

CH = 9m- 2m m?m-l- ( 2 - 53 ) 
We will return to this general strategy in a later paper and will be content to 
illustrate it with the examples in the next section. 

3. One parameter deformation of the Chebyshev polynomials 

The two dimensional Chebyshev polynomials corresponding to the measure 

\/l — x 2 \/l — y 2 dxdy 

are parametrized in terms of the Sij from Theorem l2.3l as follows so,o = 1, S2n-i.o = 

0, S2n,0 — f' S 0,2m-1 = 0, S() : 2m = \, S2n-l,2m-l = 0, S2n-l,2m = 0, S2n,2m-1 = 0, 
S2n,2m = \ for U, TO <E N. 

In order to have a nontrivial two dimensional polynomials (i.e. the measure is not 
just a product of two one dimensional measures) it was shown in [5] that %ij =/= 
for at least one pair of indexes (i,j) £ N 2 . One way to construct such polynomials 
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with relatively simple recurrence coefficients is as follows. We take 3Ci,i = [si,i] 
to be a nonzero matrix, i.e. pick si,i 7^ 0. Following the algorithm in Section 6 
of [3j , we pick at each level (n, m) with 71 j 771 M parameters S2n— 1,2m— 1 — ^1,1; 

s 2n-l,2m = 0; S2«,2m-1 = 0, S2n,2ro = 2' 

This specific choice leads to a one parameter deformation of the Chebyshev 
polynomials with recurrence coefficients given by the following formulas 












: 


si,i 



7 1 = _ 

ro,m 2 



f 
f 1 



and m is a zero matrix. For A n;m and B rhm we obtain 

"f 



1 



1 



R - 
-t>0,m — ~ ^ 



f 

1 f 



f f 
f 



and for n > 2 we have 



f 

f 



%-Lm — 0, 



(3.1) 



i.e. beyond n = 1 the matrices ^4„ im and B n ^ m reach their asymptotic values. 
Similar formulas hold if we reverse the roles of x and y, i.e. we have a measure 
symmetric in x and y. 

Equation (|3.ip shows that the right-hand side of formula (|2.23p will vanish for 
n > 2, which gives that 

K,m( Z ) = **,m(*) fOT " > 1. 

and f+(z) = j^^l m by (|2.28p . This gives the matrix measure dM(x) via formula 
(I2~35)l . 

The entries of the Hankel matrix dM m+1 (x) allow us to compute the moments 
Hj{x)dx = J y J d(j, y (x,y), thus the coefficients in the three term recurrence relation 

For the first few values one obtains that 



'1,1 



Sl.lX 



and 



(3.2) 
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This suggest one possible way to compute the measure d/j,(x,y), by first proving 
that (|3.2p holds and then using one dimensional scattering theory in the variable 
y with a; as a parameter. This would give that w m ^pm( w ) i s independent of m for 
to > 1, where is the function defined by (|2.53[) . i.e. 



Sf jU) — 2si ixw + 1 



1 - s 2 



Applying the one dimensional scattering techniques, we can show that the measure 
is given by the following formula 



dfj,(x,y) = — ^ — V y dxdy 

4 (1 - sDv 7 ! - ^v 7 ! - F 



(3.3) 



7T 2 4s 2 x (a; 2 + y 2 ) - 4s M (l + s\ Jxy + (1 - s\ J 



;dxdy. 



In our particular case, we can give a simpler proof using the connection of the 
polynomials above with the Chebyshev polynomials. This is the content of the 
next theorem. 

Theorem 3.1. The polynomials ^? n ,m{x,y) constructed with the parameters Sij 
given by 

s 0,0 = 1) s 2n-l,0 = 0, S2 n ,o = ^' s 0,2ro-l = 0, s 0,2m = 
S2n— l,2m-l = s l,l, *2n— 1,2m = 0, S2n,2m-l = 0, S2n,2m = 

for all to, n € N are orthonormal with respect to the measure dfx(x, y) defined by 



Proof. We have 

lPn,m — Pn,rn{x) 

Thus the orthogonality relation 



(3.4) 



i(x,y)Pl (x,y)dn{x,y) = 8 nyk I, 



m+l 



(3.5) 



is equivalent to 



\ 



[I y ... y™]o> 



where 



dy/i 



J 



2 V / r^ 2 



Pk,m {x) 



dfi 



(3.6) 
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is a positive measure in the variable y. From (|2.35p we know that 
f l 0\/l — r 2 

p^M^Ti^^Ti^r'PkMx)— dx = 5 n , k I m+1 . (3.7) 



If we denote 



(3. 



then to prove the orthogonality relation (|3.5j) , it is enough to show that 

0i,m(2, y)4>i, m (z, y)^dyfi = I m +i- (3.9) 



The proof of the last equality can be easily established by using two facts. The first 
one is to notice that 



i.e. for the measure given in 



o2 _2 



2si Azy + 1 



we have 



2 v/i^y 5 



(3.10) 



The second step is to connect (j>i, m (z, y) to the Chebyshev polynomials of the second 
kind. More precisely, we will show that 

4>i,i(z,y)Ui{y) if i<m 

\U m (y) - sx,iU m -i(y) i£i = m. 

Notice that l/</>? i(z, y) is essentially a generating function for the polynomials 
U m {y). This combined with elementary properties of the Chebyshev polynomials 
leads to 



<P\,m(^y) 



(3.12) 



^m(z.!/) = e>i,i4>i,i{z,y)^2(za 1 ,i) i U tlH . j (y). 

j=o 

Equation ([33)) follows immediately from (|3~TU)) - ([3~T^|) . 

Thus it remains to establish (|3.1ip . From the recurrence coefficients, it is clear 
that 

pd,m( x i v) = Poj( x , v) = u j{y) for ever y i e N o- 

From (|2.1ip with n = 1 it follows that 

Pimfo 2/) = Pl.m-lfo 2/) for < J < 772 - 2. 

This implies that 

^l.mO*) y) = ^i,m-i( z > 2/) for < j < m - 2. 
Since m = —% ntm A t n _ l m relation (|2.12p for n = 1 gives 



(3.13) 
(3.14) 
(3.15) 



A A + Pl + rl = 2 OTl,m-l = 2 yPl,rn for J < m ~ 2 



i _ 2 n m.-2 , m 



%2>™m-i - si,iZ7 m _i(y), 



(3.16a) 
(3.16b) 
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where in the first equality we used (|3.14|) and in the second (13.13p . Using the 
definition of 4> n , m , (I3.i3|) and (|3.f6a,p it follows that 



y l,m 



= tyftl m for 3 < m - 2. 



The last equation combined with the definition of 4>i,m leads by induction to §\ = 
4>\ \Uj(y) for j < m — 1 thus proving the first part in equation (|3.f f p . 

Finally, beginning with the definition of </>"\„ and then using equation (|3.f 6bj) 
to eliminate p™m yields 



2^1^-1 ~ *l,lE/m-l(lO - \/l - «? iP^ 2 



zU m (y). 



Now substituting in the definition for 0™m-i' using the recurrence formula for 
2yU m -i(y) and the definition for <^™~ 2 gives 



ro-2 



from which it follows 



2ysj 



U m -2(y) - si A U m -i(y), 



which gives the second part in p. lip . 



□ 



Remark 3.2. We have computed also the first few coefficients for the recurrence 
relations (|1.1|) - (|1.2|) in the total degree ordering. The computations suggest the 
following formulas 



A — i 



Sl,l 



s 2 
'1,1 



(3.17) 



A — 
■"■V,n 



I o 



(3.18) 



and B x ^ n , B y ^„ are (n + 1) x (n + 1) zero matrices. 



4. TWO PARAMETER DEFORMATION OF THE CHEBYSHEV POLYNOMIALS 

Another deformation which leads to simple recurrence coefficients can be ob- 
tained by making the same choice of parameters as before except at levels (1,0) 
and (0, 1). At level (1,0) we leave si : o free and we choose S2.0 — \- At level (0, 1) 
we put so,i = si,osi,i and so, 2 = \- Using the same algorithm, we have this time 
two free parameters s\ t i and si.o and the recurrence coefficients take the form (we 
denote them by A', B' , etc. in order to distinguish them from the one parameter 
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deformation) 



A' 



1 - s 2 



1 

—281,081,1 1 



-2si j0 si,i 1 

-2si si 



A 'n,m = |-Wl for n > 2 , 



R' — 



2si, si,i 

si,i 2si, (l - s?i) si 



Si,l 251,0(1-81,1) Sl,i 

si,i 2si, (l - s 2 ^) 



s l,0«i 1 



Sl.oSl.l 



Sl.OSi 







B' = for n > 2, 



and the coefficients 3C„ !m and J* m are the same as in in the one parameter defor- 
mation. 

Theorem 4.1. The polynomials V n m (x,y) corresponding to the parameters Sj,j 



given by s ,o = 1, |si,o| > h, s 2 ,o 



S2n-1,0 



0, s 2 nA = h, for n > 2, 



So,i = 81,0814, s 0: 2 = \, s 0: 2m-i = 0, s ,2m = \ for m > 2, and 

S2n-l,2m-l = s l,l, S2 n -l,2m = 0, S2n,2m-1 = 0, S2n,2m — Tl, 771 G N 



are orthonormal with respect to the measure 

2z 



dfi'(x,y) 



■K 2 {X ~ xj ° 



(x, y)yl — x 2 \J\ — y 2 dxdy 



— — — Mo(#, y)S{x - x )s/T 



(4.1) 



y 2 dxdy, 



where zq — l/(2si,o) an d 
Ho{x,y) = 



1,1) 



2 \2 - 



(4.2) 



4s 2 1 (x 2 + y 2 ) - 4si,i(l + fl? iX )xy + (1 - s 2 ^) 

Proof. The shortest way to derive the orthogonality measure is to use Theoremf 
Indeed, let L m be the second order difference operator corresponding to the second 
order difference equation (|2. 10[) for the polynomials P n ,m(x,y) defined in Theo- 
rem EH1 i.e. (£ m (/))„ = A n +i, m fn+l + B n>m f n + A^^fn-i. Similarly, we de- 
note by L'm the operator corresponding to the polynomials f' n m (x,y) defined in 
Theorem 14.11 Then one can check that L m and £/ m are related by a Darboux 
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transformation if we take xo — ^(zq + 1/zq) — (4s 2 + l)/(4si i o) and the operators 
IP and Q are defined as follows 

nf)n = \f n+ i - -r^—U for n > 1 and 3>(/) = Ai, ro /i ~ -^—(a^fo 

Q(/)n = In - 4s l!0 A^ m /„_i for n > 1 and Q(/) = a m f 0l 
where a m is the (to + 1) x (to + 1) matrix given by 

1 

-2si,oSi,i 1 

-2si Si l 1 



-2sinSi 1 1 



From this relation it also follows that 



(4.3) 

where P„ jrn are the polynomials defined in Theorem 1 2. 41 Writing the orthogonality 
relations for P n , m (x), using (|4.3|) and replacing P^ m (x)[l y . . .y" 1 ] 1 by F' nm (x, y) 
we obtained the desired orthogonality. □ 

Remark 4.2. Finally we list the recurrence coefficients in the total degree ordering. 
The matrices A x<n and A y ^ n are the same as the matrices for the polynomials in 
Section 3, i.e. they are given by (|3.17p and (|3.18|) . respectively. For the B x n we 
have 



Bx.o = l s i,o]> B Xt i = si.o 
and for n > 2, B XjTl is the block matrix 

Bx.n = 



1 - S 2 



-si,i A /l - s{ 



B m ,i 




'1,1 



The matrices B y ^ n are identically equal to zero for n > 1 and B y fi — [si^si i]- 
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